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The asymmetric Hubbard model is used in investigating the lattice gas of the 
, moving particles of two types. The model is considered within the dynamical 

£h ' mean-field method. The effective single-site problem is formulated in terms of 

the auxiliary Fermi-field. To solve the problem an approximate analytical method 
based on the irreducible Green's function technique is used. This approach is 
tested on the Falicov-Kimball limit (when the mobility of ions of either type 
is infinitesimally small) of the infinite-?/ case of the model considered. The 
dependence of chemical potentials on concentration is calculated using the one- 
1 particle Green's functions, and different approximations are compared with the 

J"""-. , exact results obtained thermodynamically. The densities of states of localized 

' particles are obtained for different temperatures and particle concentrations. 

The phase transitions are investigated for the case of the Falicov-Kimball limit 
in different thermodynamic regimes. 
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a 

o '■ 1- Introduction 



PACS: 71.10.Fd, 05.30. Fk, 05.70.Fh 



Strongly correlated electron systems have been a subject of growing interest in 
recent years. These systems became especially interesting after the discovery of high- 
T c superconductivity. In general, the methods used for their description are based on 
the Hubbard model and its generalizations taking into account strong short-range 
correlations of particles. Similar models are considered in investigating the ionic 
conductivity in crystalline materials. Fermi lattice gas models can be mentioned 
among them. The asymmetric Hubbard model [1] arises by extending the models 
to the case of the systems with moving ions of two types. The transfer parameters 
and chemical potentials are different for the ions of different nature. As some spe- 
cial cases, the asymmetric Hubbard model includes the Falicov-Kimball model and 
the standard Hubbard model. Thereby, it should be mentioned that the asymmet- 
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ric Hubbard model has been recently obtained as a generalization of the electron 
Falicov-Kimball model which describes the interactions between itinerant (^-electrons 
and localized /-electrons by including /-/ hopping [2]. The model was used in de- 
scribing the effects related to orbital ordering, as well as to Bose-condensation of 
electron-hole pairs (which includes a spontaneous ferroelectric polarization) in oxide 
compounds. 

Despite the relative simplicity of these models, the theory of electron spectrum 
and thermodynamic properties of such systems is far from its final completion. In 
recent years the essential achievements of the theory of strongly correlated electron 
systems have been connected with the development of the dynamical mean-field 
theory (DMFT) [3]. This method is exact in the limit of infinite space dimension 
(d = oo) and is based on the fact that the irreducible (according to Larkin) part of 
the one-particle Green's function is single-sited (at appropriate scaling of transfer 
parameter) . 

The central point in this method is the formulation and the solution of the 
auxiliary single-site problem. In this case a separated lattice site is considered as 
placed in some effective environment. Since the processes of the particle hopping 
from the site and returning into the site are taken into account, the mean field 
acting on the particle states of the site is of a dynamic nature. The field is described 
by the coherent potential J a {uj) that should be determined in a self-consistent way 
from the condition that the same total irreducible (according to Larkin) part defines 
Green's functions both in the atomic limit and for the lattice [4-6]: 



Here a is the type index (or a spin index for electron systems); there is a possibility 
of the transfer parameter t£ can depend on a particle type. Summation over the 
wave vector can be changed by integration with the density of states. There are 
cases that are usually investigated in the d = oo limit: the hypercubic lattice which 
is the generalization of the cubic d = 3 lattice, and the Bethe lattice which is the 
thermodynamic limit of the Cayley tree when the number of the nearest neighbours 
tends to infinity. The density of states is Gaussian for a hypercubic lattice [3], while 
on the Bethe lattice the density of states is semielliptic [4]. 

The single-site problem can be solved analytically only in some simple cases. 
In general, the application of numerical methods turns out to be necessary. In the 
case of the Falicov-Kimball model it is possible to investigate the thermodynamic 
properties of the system by analytical calculation of the grand canonical potential 
[7,8]. Also, the analytical expressions for the itinerant and for the localized electron 
Green's functions have been obtained [9,10]. However, the topology of the localized 
electron band has not been completely investigated yet. 
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Recently, an approximate analytical approach to the solving of a single-site prob- 
lem has been developed [11,12]. This approach is based on the irreducible Green's 
function technique with projecting on the Hubbard basis of Fermi-operators. A sys- 
tem of DMFT equations was obtained in the approximation which is a generalization 
of the Hubbard-III approximation combined with a self-consistent renormalization 
of the local electron levels (here this approach will be called GH3). It has been shown 
that the proposed approach includes, as special number of known approxi- 

mations (Hubbard-III, alloy-analogy (AA) and modified alloy-analogy (MAA)). 

In the present work this approach is formulated for the asymmetric Hubbard 
model. Its applicability within different approximations is tested on the infinite- 
U spinless Falicov-Kimball model. The densities of states of the moving and the 
localized particles are obtained for the Bethe lattice at different particle concen- 
trations and temperatures. Dependence of the chemical potentials on concentration 
is calculated using the Green's functions. The results of different approximations 
are compared with the exact results obtained thermodynamically. Some thermody- 
namic properties (such as phase transitions and phase separations) of the system in 
different thermodynamic regimes are investigated as well. 

2. Model 

We consider the case when the lattice gas model is used in describing the systems 
with moving particles (ions or electrons) of two types. A single site can be in three 
possible states: (i) no particles, (ii) one particle of the first (A) type, (iii) one particle 
of the second (B) type. Motion of the particles can be described by the creation and 
by the annihilation operators and by the transfer parameters dependent on the 
particle type. Since the site can be occupied by only one particle of either type, the 
operators are not of the Fermi type. It is possible to investigate the system using 
the Fermi operators by formal adding the fourth state of a site occupied by two (A 
and B) particles simultaneously. In this case, the asymmetric Hubbard model arises. 
The Hamiltonian is 

# = £#< + £ t*af+af + £ 4a?+af, (4) 

* (ij) (ij) 

Hi = -HAuf - fx B nf + Unfnf, (5) 

with chemical potentials dependent on the particle type. It can be easily seen that 
Hamiltonian (4) corresponds to the Hubbard model in an external magnetic field 
with the spin-dependent electron transfer. The U term is the single-site interaction 
energy between the particles and it should tend to infinity when we wish to return 
to the initial lattice gas model. 

In the limit of large dimensions, the problem is investigated using the DMFT 
approach. The single-site problem is formulated in terms of the auxiliary Fermi-field 
[11]. Let us write the effective single-site Hamiltonian in the Hubbard operators 
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representation 

H cS = -J2 [^( Xaa + x22 )] + UX22 

+ [(X a ° + C* 2 *)^ + + CX" 2 )] + H 6 , (6) 

<J 

where the notations are used: a = B, ( = + for a = A; a = A, ( = — for a = B. 
The basis of single-site states \n A ,n B ) is 

|0) = |0,0), \A) = \1,0), ( 
|2) = |1,1), |S> = |0,1>. [n 

The Green's function G^\u) = {{a^a^))^ can be written in the following form: 

G (a) = ((x°-|X CT0 )) w + C{{X^\X 2 ")) W + a(X* 2 \X°°)) w + ((X^X 2 ^ . (8) 

The auxiliary Fermi-field describes the environment of the selected site and is for- 
mally characterized by the Hamiltonian H^. The explicit form of this Hamiltonian 
is unknown. However, to calculate the Green's function G^\u), the averaging over 
the £, £ + operators is done with the help of the function 

QM = «^IC»? e) - (9) 

The relation 

2nV 2 Q a (Lu) = J a (cu) (10) 

takes place in this case. Unlike the standard Hubbard model, the coherent potential 
J ry is dependent on the type of the particles (J a ^ Jb for tfj ^ tfj). 

3. Green's functions for the effective single-site problem 

Equations for functions (8) are written using the equations of motion for X- 
operators. According to the method developed in [13,14], let us separate (in the 
Green's functions of a higher order) the irreducible parts expressing the derivatives 
i dX °J a2) dt = [X 0<t(ct2 ), H c g] as the sums of regular (projected on the subspace formed 
by operators X 0a and X a2 ) and irregular parts 

[X 0ct , H cS ] = -fi a X 0a + al u X 0a + a° 2 a X 52 + Z 0a , 

[X* 2 , H cS ] = (U - i2 a )X° 2 + af 2 X 0a + a? X° 2 + Z° 2 . (11) 

Operators Z 0a and Z u2 are defined as orthogonal to the operators from the basic 
subspace and we come to the expressions 



Z Qu = V{X 00 + X<">)& + VX™£ 9 + c^x 02 ^, 

Z° 2 = V{X 22 + X^% - (VX^ - VX^H, (12) 
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where 



( X oo + X ™)£ a = (X 00 + X™)t a , 
(X 22 + = (X 22 + X™)^, 



X aa^_ = X aa^_ _ J_fe X oO )X 0* _ J_ ( X 2 °^)X* 2 , 

l / , ^^-O 1_ \ -m- rC\^ l 



X 02£ = X 02^ _ i x <ng )x *> _ { gx<» )x n (13) 

and A pq = (X™ + Xii); Aoa = 1 _ n . j A 2 _ = n _. 
In this case 

a°* = -Caf 2 = 

<p < , = foX°°)+C{X« l &). (14) 

Using this procedure by differentiating both with respect to the left and to the 
right time arguments, we come to the relations between the components of the 
Green's function G^ and scattering matrix P c . In a matrix representation, we have 

G a = Gq + GqPvGq, (15) 

where 

/((X 0ct |X ct0 )) {{X°°\X 2 °)) 
\((X° 2 \X C 

and nonperturbed Green's function Gq is 
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where 



D a = (LU - a a )(u - b a ) - -T—l <Pa: (18) 

V V 

-H„ + "j— <Ar, b a = U - fjL a + -r—<Pa- (19) 



The scattering matrix 



p - 27r ( A ^ \ ( ((Z°°\Z° )) {(Z»°\Z 2 -°))\ ( A£ ] 

F ° ~ i7r I o A£ ) I «z«|Z"°» «z«|z»» J a 2 1 ' UJ) 



2(7 



being expressed in terms of irreducible Green's functions contains the scattering 
corrections of the second and the higher orders in powers of V . The separation 
of the irreducible parts in P a enables us to obtain the mass operator M a and the 
single-site Green's function expressed as a solution of the Dyson equation 

G ff = (l- GlM^Gl. (21) 
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We will restrict ourselves hereafter to the simple approximation in calculating 
the mass operator P a , taking into account the scattering processes of the second 
order in V. In this case 

M a = Pj°>, (22) 

where the irreducible Green's functions are calculated without allowance for corre- 
lation between electron transition on the given site and environment. It corresponds 
to the procedure of different-time decoupling [15], which means in our case an in- 
dependent averaging of the products of X and £ operators. Let us illustrate this 
approximation by some examples. 

1. The Green's function (((X 00 + X™)£ CT |£+(X 00 + X™))) w = h(u). 

Using the spectral theorem and performing the different-time decoupling: 

(C(t)(x 00 + x™) t (x 00 + x aa )^) h » ((^ 00 + X™) t (X 00 + X™)) 

(23) 

we obtain 

h(u) = A 0a ((C \£)) u = -^JM- (24) 
Here, X-correlators are in a zero approximation 

((X 00 + X™) t (X 00 + X™)) « ((X 00 + X™) 2 ) = A 0a . (25) 



2. The Green's function {{X s * £ 9 \£tX™)) u = I 2 {u). 

The time correlation function (££(t)X c " T {t)X aa £ sd . ) ir is decoupled as 

i&wwt, r « (x^wK&Mt, ). (26) 

In the zero approximation 

(X™(t)X*°) = exp [i(^ - ^)t](X™). (27) 
Using these expressions we obtain 

r 2rr ^ (X<">) + (X 99 ) , . (X<">)-(X™) 



2irV I 2 (UJ) = Jaiyj + He - fia) ~ 



2 uy ^ u ^ 4tt 

"+00 

X 



/+oo j / o I 
^ tanh^-{-2ImJ 5 (cu , + i£)}. (28) 
U - U' - flu + fir 2 

Let us notice that in the case with ua = tib, = Hb, which corresponds to 
the simple Hubbard model in the absence of an external magnetic field, 
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3. The Green's function ({X 02 &\^ X 20 ))^ = I 3 (uj). 
In this case 

1 3{ UJ ) - ^ Iw//uH-/*<H-^-eH ^ 

_ ,-,- + , a + ,„- t / ta " h V'- 2Im «gfc»-'- < 30 > 

Using the above results, the mass operator can be obtained in an explicit form. 
According to (22), it is given by the expression (20) where 

{(Z 0ff \Z a0 )) u = A 0( ,Mu)-R ( ,(u), 
((Z^IZ 2 *))^ = A 2 - a J a (uj) - R a (u), 

((Z°*\tf*)) u = ((Z* 2 \Z° ))^ = R a {u), (31) 

and 

(X <Tcr ) + (X°°) 

Ra{u) = ~- '-^ L J S (u + Ha) 

+ - - ^ - / {-2ImJ ff (u' + ie }tanh — 

+ (x° ) + (X a ) MU _ tia _ ti9 _ u) 

+ ~ hr^ " / ^ ; tt{- Im M-u' - is)} tanh H 



2tt 7-oc ^ - ^' + /% + /V - U " " 2 

(32) 

4. Set of DMFT equations in the infinite-^/ limit 

The following consideration will be performed in the case of U — > +oo which 
excludes simultaneous occupation by two particles of A and B types of the same 
site, when the model is used in describing the lattice gas of the particles of two 
types. Then functions R a for a = A or B are as follows 

R a {uj) = - J B {uJ + Ha- Hb) 

. n A -n B f + °° du' (3u' , . 

H ; / ; tanh— — { - 21m J B (uo + ie)\, (33) 



Rb(oo) 



+ 



n A + n B 
2 

n B - n A 
An 



J A (u + jj, B - jJ, A ) 



f3u' 



tanh-— { -2ImJ A (cc/ + ie)}, (34) 

u - uj' - ji A + hb 2 
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The single-site Green's functions (8) are obtained using relations (21), (22) and (31) 
in the U — > +oo limit, and respectively they equal to: 



G ( a\^) = * HB ^TT, (35) 



G b(") = — / HA T( ^ (36) 

Parameter </? CT , which is expressed by the average values of the products of X 
and £ operators, is a functional of the potential J a {uj). According to the spectral 
theorem 

^ JTT [V((£ (J |X CT0 )) w+i£ -\/((£ (T |X^))_ i£ ]. (37) 

The Green's functions are found using linearized equations of motion and neglecting 
the irreducible parts. 
In this case 

1 /•+" da, f (l-n A )J B (u>) ■) 

Vv *=-2iL ^Til- 2h , +to -^ U,- (38) 



The coherent potential J^u) is self-consistently determined from the equations 
'l)-(3) by eliminating the total irreducible part. Integration with the semielliptic 



A^) = ( 4 °) 



density of states is done. In this case, we have simple expressions 

GfM = 4 ~^f, (41) 

= (42) 

The set of equations becomes a closed one by adding expressions for average particle 
concentrations obtained with the help of the imaginary parts of the Green's functions 
(i.e. interacting densities of states): 

1 /■ +0O J -21mG^\co + ie) 

tia = — — / duj -, (43) 

^ Loo e^ + 1 ' K J 

1 -2ImG ( R a) (w + ie) , . 

2vr/ e^ + 1 1 ; 
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It has been shown in the case of the standard Hubbard model [11] that such a 
set of equations corresponds to GH3 approximation, and includes, as some special 
cases, a number of known approximations. The case with R a = 0, <p a = corresponds 
to the alloy-analogy (AA) approximation [16]. The modified alloy-analogy (MAA) 
approximation is obtained by taking into account the renormalization of the local 
electron levels tp a [16,17]. For R a ^ 0, ip a = the system of equations corresponds 
to the extension of the Hubbard-III approximation by the inclusion of the integral 
terms responsible for scattering. 



5. Lattice gas thermodynamics in the Falicov-Kimball limit 

5.1. Exact solution of DMFT problem 

In the limiting case of U — > +oo, W B — > (tfj — > 0), when mobility of the B type 
ions is very small (this case corresponds to the spinless Falicov-Kimball model [18] 
with the infinite on-site repulsion), the effective single-site problem can be solved ex- 
actly. Thermodynamics of the model, equilibrium states and phase transitions were 
investigated in a series of papers ([7,8,19], see [20] as well). There were considered 
thermodynamic regimes with the constant n = ua + »b, £ = Hb — Ha values [8] and 
with the fixed relative filling % = const, ha = const • (1 — [7]. However, other 
thermodynamic regimes should be also considered: (i) constant chemical potentials 
Ha, Hb', (h) the mixed regimes ha = const, n B = const or \i B = const, ua = const. 
Let us note that in the case of the pseudospin-electron model (PEM) without tun- 
nelling, the dynamics of pseudospins being analogous to the Falicov-Kimball model, 
the regime of hb = const corresponds to the fixation of longitudinal asymmetry 
field h acting on pseudospins. The PEM thermodynamics has been a subject of in- 
vestigations within DMFT [21] as well as in the framework of approximate methods 
(generalized random phase approximation (GRPA), [22-24]). 

An exact solution of the single-site problem is given by the simple expression: 

G { l\u) = , v (45) 

It should be mentioned that this formula can be obtained by the above described 
approach and corresponds to the AA approximation. The final expression for the 
Green's function is obtained by solving the quadratic equation formed by (41), (45): 



G ( 2'(uJ±ie 



uj _i_ 

2 <u + ha) + w^Vi^ + Va) 2 - W 2 A (l - n B ), uj + ha< -W Ay /l - n B ; 



Wa , 

^(oo + ha) T i4i^ W a( 1 - "flj - + Ha) 2 , \u + Ha\ < W A ^/l - n B ; (46) 

2 

w\y~ ' Wa 



{uj + fi A ) - T^y/(u + Ha) 2 - W\{1 - n B ), uj + fi A > W A y/l - n B . 



Here, the phase of the square root is chosen so that an imaginary part of the Green's 
function has the correct sign as well as by using the properties which follow from 
the spectral theorem. 
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Figure 1. (fiA, /Ufi)-phase diagrams. The 
parameter values Wa = 1, = 0, the 
dashed line - T = 0.01, the solid line - 
T = 0.04, the bold line - T = 0.06. 



Figure 2. The part of (/xa, /U_B)-phase 
diagrams at different temperatures; the 
appearance of the phase transitions for 
HB < (Wa = 1, W B = 0) is shown. 



Using this Green's function the first equation for the particle concentrations is 
obtained from (43) 



20-ng) f +1 VT^ .... 

n A = 1 / drr — — — — (47) 



Here, in the standard scheme, the second equation is obtained thermodynamically 
by differentiating the grand canonical potential [12] 



where in our case 



n - r dx v 7 ^ 1 (49) 

^~ ^ 7T e -P(xW A VT^-»A)) + 1 arctan x • l4yj 

Let us notice that parameter Q does not depend on the sign of chemical potential 
of the moving particles, i.e., it does not depend on the interchange ha —^a- 

Separating the solutions of a set of equations (47)-(49) that correspond to abso- 
lute minima of the grand canonical potential, we will investigate equilibrium states in 
the above-mentioned thermodynamic regimes. There are phase transitions between 
homogeneous phases with different particle concentrations in the regime of constant 
chemical potentials (/ia = const, /ib = const). Phase separation phenomena take 
place in the regimes where either of the particle concentrations is constant. 

1. ha = const, hb = const. In figure 1 phase diagrams are shown in (/ia, I^b) 
coordinates. All energy quantities are given in the units of a half-bandwidth 
(Wa = 1). The phase transition curves terminate at critical points. They 
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(a) (b) (c) 

Figure 3. The phase separation diagrams for Wa = 1, Wb = 0. (a) - (ub,T)- 
phase diagrams for the different chemical potentials of moving particles, (b), (c) 
- (ha, T)-phase diagrams for /xb = —0.006 and hb = 0.6 respectively. 



become shorter (symmetrically with respect to the value \ia = 0) with the 
temperature growth and vanish at critical temperature. Numerical calculations 
give the following values of critical parameters 

T c ^ 0.0601W^, a4 = 0, fi% fa Q.229W a . (50) 

At the zero temperature, phase transitions are within the —Wa < Ha < Wa, 
< us < Wa ranges of chemical potentials. However, figure 2 shows that there 
is a possibility of phase transitions with \xb < when temperature increases. 

Let us note that diagrams analogous to those in figure 1 can be obtained in 
DMFT for PEM (where ha corresponds to chemical potential fi of electrons, 
and hb corresponds to the asymmetry field h). The phase coexistence curve 
on (/i, h) diagram obtained for finite values of the coupling constant in [21] 
consists of two parts corresponding to the chemical potential being in the 
vicinity or inside of either electron subband (the band splitting in a spectrum 
is caused by interaction with pseudospins) . There is a direct correspondence 
between the lower part of the phase diagram for PEM and the diagrams in 
figure 2. 

2. fiA = const, tib = const (figure 3(a)). In this case, a phase separation takes 
place. A topology of the phase diagrams (n B ,T) depends only on the absolute 
value of the chemical potential of the moving particles. The highest value of 
critical temperature is reached at /ia — and is given by the expression (50). 
The critical temperature decreases when chemical potential of the moving 
particles approaches either of the band edges. At the zero temperature for 
—Wa < Ha < Wa, the system is separated into two phases: = and 

(2) 

n B = 1. The regime of relative half- filling (% = const, ua — (1 — ^s)/2) has 
been considered in [7]. This case corresponds to the % = const regime with 
Ha = 0, so the diagram which is equivalent to figure 3(a) for /ia — has been 
previously obtained. 
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Figure 4. The ub dependence of fis in different approximations is compared with 
the exact result obtained thermodynamically. The parameter values: Wa = 1, 
W B = 0. 1 - exact result; 2 - AA; 3 - MAA; 4 - H3; 5 - GH3. 



3. ua = const, Ub = const (figure 3(b,c)). There is a homogeneous state at 
zero temperature, but a segregated state appears at higher temperatures for 
—0.009Wa < Ub < (figure 3(b)). The phase separation region decreases 
continuously and there is a point at temperature T pa 0.01 Wa when chemical 
potential of localized particles approaches some critical value \ib ~ —0.009 Wa- 
The possibility of such a situation is illustrated by the behaviour of the left 
end of the phase transition curve in the vicinity of /ib = 0, /ia = —Wa 
values (figure 2). The thermodynamic regime with e° = ha — /Ab = const, 
tia ~\~ f^B = const has been investigated in [8], and phase diagrams similar to 
that in figure 3(b) have been obtained. It can be seen that both ends of the 
phase transition curve (ha,Hb) have the similar temperature behaviour (with 
respect to jj, B = and to /i# — /ia = lines). 

At /i B > 0, the topological structure of phase diagram (n A) T) is similar to the 
one in figure 3(c). A phase separated state always exists for low temperatures 
with non-integer values of ha in both phases at T = 0. Diagrams of this type 
have been obtained for PEM in the strong coupling case (that corresponds to 
U > W A ) using a GRPA approach [23]. 



5.2. Approximate analytical approach 

Let us consider the results obtained at the Falicov-Kimball limit of the considered 
model with U — > +oo in the case when the above described approximate approach 
is used instead of the exact one (48). In the Falicov-Kimball limit, Jb = and it 
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leads to (p A = 0, R A = 0. According to equation (44) the argument in the expression 
for R b {oj) (34) should be replaced by u + ie. Substituting the Green's function given 
by (46) into this expression we obtain 

r> / i • \ n A + n B 

R B (u + le) = J A {u + is + p B - n A ) 

..n B -n A . (3(u - n A + p B ) 
+ l lm J a {uj + 15 + fi B - p A J tanh 

n B -n A 2 f 1 dxVl - x 2 tanh § {xW A y/l - n B - p A ) 

+ — ~ A W A (l-n B )Vp — , 51 

47T }_ x UJ + fl B - xW A y/l - n B 

where the principle value of the integral is taken. The similar procedure for equation 
(39) leads to 



VfB = - ^ (1 -" b)3/2 Vp /' ixVY ^ [l + ^»^-«) 
27r y_! X L 



_1 . (52) 



The R b (uj) function provides in this approximation a frequency dispersion of the 
Green's function of localized particles and leads to the broadening of the correspond- 
ing density of states into a band, (which is shifted with respect to the initial level 
on the — ^ <Pb value). 

The dependence of the chemical potential on the concentration of localized par- 
ticles is calculated using expressions (44), (47), (51), (52). In figure 4 the approxi- 
mations such as AA, MAA, GH3 and H3 (in the last case the integral terms that 
describe the processes of scattering by boson excitations are neglected in the ex- 
pressions for R a ) are compared with the exact results obtained thermodynamically. 
The simple AA approximation cannot describe phase transitions because it gives a 
monotonically increasing dependence \i B = fi B (n B ). MAA and AA approximations 
can describe the thermodynamics of the system at high temperatures or at a neg- 
ative value of /jL A (small concentration of moving particles). The H3 approximation 
gives better results for high temperatures but not for the low ones. The best results 
are given by the GH3 approximation for high temperatures, as well as for low tem- 
peratures in the cases of a nearly full or empty band of the moving particles (and 
high values of n B ). 

In figure 5 the case of a nearly full band of the moving particles is shown. The 
fJ>B = /Ub(^_b) dependence is illustrated together with the corresponding curves for 
density of states of localized particles in H3 and GH3 approximations: 

p B (iu) = --lmG ( B \uj + ie). (53) 

7T 

One can see that this band can have a complicated structure or even be split in 
the H3 approximation, but it does not correspond to the satisfactory /i B (U B ) de- 
pendence. In general, there is a considerable difference between densities of states 
p B {oj) given by the H3 and GH3 approximations, but the band edges are determined 
properly even in the H3 approximation. 
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Figure 5. /ig as a function of ub in different approximations (1 - exact result; 2 
- AA; 3 - MAA; 4 - H3; 5 - GH3) and corresponding densities of states ps of 
localized particles (1 - GH3; 2 - H3). 



In figure 6 the densities of states of localized particles are presented for different 
temperatures and particle concentrations. In the cases corresponding to the above- 
mentioned criteria ((i) high T; (ii) low T and chemical potential ha close to the 
edges of the band of the moving particles), the obtained densities of states ps can 
be considered as close to the real ones due to a satisfactory description of the /^(ne) 
dependences. 

Calculation of densities of states of localized particles in the Falicov-Kimball 
model had been done earlier within DMFT in [9,10]. The algorithm which leads to 
exact results is rather complicated. The plots obtained in [9,10] correspond to finite 
values of U and to the symmetrical case of half-filling (the hypercubic lattice with 
Gaussian density of states was considered). The band of localized particles is split 
and consists of two subbands for U larger than some critical value (U > 5 Wa)', 
for U — > +oo, the only lower subband remains. Our results correspond by the 
general shape of the curves to the above-mentioned results in the cases illustrated 
in figure 5 for T = 0.02 0.2; ha = 1-2; % = 0.9 and in figure 6 for T = 0.2; 
Ha = 0; ub = 0.9 within the GH3 approximation (in these cases, there is the best 
agreement with the exact results for p>b{jib) dependences). Quantitative comparison 
of densities of states of localized particles will be possible after doing calculations 
with nonperturbed density of state like in [9,10] and for finite values of U . 

6. Conclusions 

Application of the approximate analytical method of solving the effective single- 
site problem within DMFT gives a basic set of equations for evaluating single-particle 
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<a+\i B co+|i B cim-h b 

Figure 6. The approximate densities of states of localized particles (1 - GH3; 
2 - H3). 



Green's functions and for investigating an energy spectrum of the asymmetric Hub- 
bard model describing the system of particles of two types with different transfer 
parameters (t A ^ t B ) on a lattice. 

The limiting case U — oo (exclusion of a double occupation of a site) and t B — > 
(the case of the inffnitesimally small mobility of particles of one of the type, when 
the model becomes the Falicov-Kimball model) is considered. 

Phase transitions are investigated in thermodynamic regimes of the fixed values 
of the following quantities: (i) ha, I^b] (ii) (J>a, ns', or (iii) ua- Phase diagrams 
describing coexistence curves of homogeneous phases with different iia and n B val- 
ues and phase separation regions are obtained. These results extending a space of 
thermodynamic parameters complement the existing information about the thermo- 
dynamics of the Falicov-Kimball model in the U — > +oo limit. 

The set of self-consistent equations, which relates the concentrations (chemical 
potentials) of particles of A and B types, coherent potentials Ja,b(v) an d shift 
constants of energy spectrum <£a,b, is solved in two cases: (i) with an application 
of the exact thermodynamic relations for the Falicov-Kimball model to derive an 
equation for concentration % of localized particles; (ii) by using the approximate 
analytical scheme in DMFT. Calculation of //b(wb) dependences being done by the 
two mentioned methods allows us to investigate an applicability of this approximate 
approach. The results closest to the exact ones are obtained within the GH3 ap- 
proximation for high temperatures, as well as for low temperatures in the cases of 
nearly full or nearly empty band of the moving particles and for high values of ub- 

Densities of states Pb(w) of localized particles are calculated using the approxi- 
mate DMFT scheme. The results obtained in the GH3 approximation can be con- 
sidered as the closest to the real ones, when a satisfactory dependence of ^(ng) is 
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achieved. Extending this approach to the cases of low temperatures and intermediate 
band filling requires an improvement of the applied scheme, namely, a more precise 
(probably self-consistent) way of taking into account the processes of scattering by 
boson excitations for calculating the Rb{oj) function which determines the energy 
spectrum of localized particles. 
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